Statics and Inertial Dynamics of a Ruck in a Rug 
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We consider the familiar problem of a bump, or ruck, in a rug. Under lateral compression, a rug 
bends out of the plane forming a ruck — a localized region in which it is no longer in contact with the 
floor. We show that when the external force that created the ruck is removed, the ruck flattens out 
unless the initial compression is greater than a critical value, which we determine. We also study the 
inertial motion of a ruck that is generated when one end of the rug is moved rapidly. We show that 
the equations of motion admit a travelling ruck solution for which a linear combination of the tension 
and kinetic energy is determined by the ruck size. We confirm these findings experimentally. We 
end by discussing the potential implications of our work for the analogous propagation of localized 
slip pulses in the sliding of two bodies in contact. 

PACS numbers: 87.10.Pq, 46.40.Cd, 91.55.Fg, 61.72.Hh 



As well as being something of an annoyance in every- 
day life, the familiar bumps (or rucks) that form in a 
rug as people walk over it has long proved to be a use- 
ful analogy in explaining a range of important physical 
phenomena. For example, several authors have used the 
motion of a ruck to illustrate how dislocations can fa- 
cilitate the relative motion of two crystalline planes of 
atoms (see pg. 117 of [lj). In this situation the ruck is 
quasistatic. However, those working with carpets every 
day know that shaking one end leads to the propagation 
of rucks, which facilitate small-scale sliding [2|] . This dy- 
namic scenario is used to explain the observation of Schal- 
lamach waves Q: waves of detachment that control the 
sliding of rubber surfaces 0, 0| ■ In much the same way 
the slip pulses observed in some earthquakes @, 01 are 
often thought to be akin to moving rucks in a rug [8( . 

Though these analogies are frequently referred to in the 
literature, we are not aware of any investigation of their 
quantitative validity nor even a study of the properties 
of rucks in rugs themselves. In this Letter, we consider 
the statics and inertial dynamics of a ruck in a rug. Our 
aim is to characterize the model system as a first step 
towards understanding the validity of these analogies. 

Recently, a great deal of attention has been focussed 
on understanding the statics and dynamics of thin elas- 
tic objects, of which a rug is a particular example. Static 
localized structures qualitatively similar to a ruck in a 
rug have been observed in the compression of thin films 
floatin g on water 0] or deadhering from a polymer sub- 
strate [10j. Generally, studies of the dynamics of thin 
elastic objects have been focused on investigating the in- 



teraction of a flexible object with a fluid [ll|, LL2J. The 
purely elastic propagation of localized disturbances along 
thin objects has received less attention, notable excep- 
tions being investigations of a wave in a whip II 311 - the 
travelling waves that form in suspension bridges [14| and 
the rolling mode in vibrated thin plates [HI]. 



We begin by considering the static problem of a two- 
dimensional sheet of thickness h, density p and bending 
stiffness B lying on a rigid, horizontal substrate. A cer- 
tain end-end displacement Al is then imposed symmet- 
rically (see fig. [I} . Since compression of the sheet itself 
is energetically expensive, we expect that the sheet will 
buckle out of the plane. Opposing this is the weight per 
unit area of the sheet, pgh, which makes it energetically 
expensive for the sheet to lose contact with the substrate 
everywhere along its length (in contrast with the conven- 
tional elastica [l6|). Instead, contact is lost only over a 
localized region — a ruck is formed. We denote the shape 
of the ruck by [x(s), w(s)] in which s is the arc-length, 
though it is more convenient to determine an intrinsic 
equation for the local inclination of the ruck, 0(s). This 
can be found by minimizing the value of the functional 

r l/2 r r n 

T = J | \B6 2 S + pgh w - Y,\Al/l- cos 6>J j ds, (1) 

in which the first two terms are the bending and gravita- 
tional energies of deformation (per unit width), respec- 
tively, and the inextensibility constraint is enforced by a 
Lagrange multiplier X — the compressive force per unit 
width in the ruck. Here I is the total arc length contained 
in the ruck and we may consider only the interval [0, 1/2] 
by symmetry. Integrating by parts and using w s = sin9, 



we find that / ' w ds = j Q (1/2— s) sin9 ds. We may 
then use the Euler-Lagrange equations to find a differen- 
tial equation for the 9(s) that extremizes the functional 
in |T]). We first note that a natural lengthscale 
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arises. This elasto-gravitational length measures the typ- 
ical sheet length required for the sheet to be deformed by 
its own weight. Introducing non-dimensional variables 
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FIG. 1: The ruck shape observed in a sheet of natural rubber 
with thickness h — 0.75 mm for different imposed end-end 
compressions AL Here the (non-dimensional) compressions 
are AL = Al/£ g = 0.3, 2.7 and 4.2. 



S = s/i g , AL = Al/i g , L = l/i g and a = T,i 2 g / B the 
differential equation for 0(S) becomes 

9 SS = -a sin 9+(L/2-S) cos 9, 



in which subscripts denote differentiation. Eqn. (l3l is 
commonly referred to as the heavy elastica equation [171 ] 
- it is the classical elastica equation [16] supplemented 
by a term representing the weight of the material. The 
differential equation ([3]) is of second order but contains an 
unknown eigenvalue a so that three boundary conditions 
are required. These are 



0(0) = S (O) = 0, 6(L/2) = 0, 



(4) 



which express that the sheet is horizontal to the sub- 
strate and torque free where it first touches the substrate 
(at S = 0) as well as the symmetry of the sheet about 
S = L/2. The boundary value problem ©-Q may read- 
ily be solved numerically using, for example, MATLAB's 
bvp4c routine fLsj]. Some example profiles are shown in 
fig. [1] superimposed upon images of experimentally ob- 
served rucks in which a fixed end-end displacement AL 
is imposed. We observe very good agreement between the 
theoretically predicted and experimentally observed ruck 
shapes though for larger values of AL there is a symme- 
try breaking instability (l9| . Although the ruck shapes 
shown in fig. [T] are qualitatively similar to those observed 
in the classical elastica [16J, a quantitative comparison 
shows that the heavy elastica shapes are flatter with a 
more localized 'bump'. This is a result of the additional 
boundary condition #s(0) = in the heavy elastica case. 

Analytical progress may be made by linearizing for 
small deformations, 9 -C 1. The resulting problem may 
be solved analytically giving 

a « 4a 2 L~ 2 , L « (768/5) 1/7 a 4 / 7 ' AL 1 ' 7 , (5) 

where a 4.49341 is the smallest positive solution of 
tanx = x. Thus for a given displacement AL we have 
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FIG. 2: (Color online) Regime diagram showing the regions 
of static friction-length (pLo), compression (AL) parameter 
space for which a ruck is observed to stick (i.e. remain) or 
slip (i.e. flatten out) when the external compressive force is 
removed. The numerically computed boundary between stick- 
ing and slipping (solid curve) compares well with the asymp- 
totic result ([6]) (dashed curve) for AL <C 1. The unit of length 



(3) is the elasto-gravitational length l g , defined in ([2]). 



the laws a w 3.441AL~ 2/7 , L w 4.844AL 1/7 and 6 w 
1.326AL 4 / 7 where S = d/l g is the dimensionless ruck 
height. These laws were previously presented as functions 
of a [2(|. Recently, energy arguments for the scalings 
alone were also given [2lj . The dependence of L on AL is 
distinct from the classical elastica where L is the system 
size and hence independent of AL. 

We next consider the question will a ruck remain once 
the compressing force that formed it is removed? For a 
stationary ruck to stick, the compressive force within the 
ruck, tr, must be smaller than the maximum static fric- 
tional force that can be exerted on the heavy sheet by 
the substrate. With a coefficient of static friction fx, we 
therefore require that a < jj, [(^o - L)/2+Nq\ in which Lq 
is the total (non-dimensional) length of the sheet and Nq 
is the normal reaction on the sheet at the point of contact 
between the sheet and substrate (i.e. at S = 0). Consid- 
ering ([3]) we see that Nq = L/2, i.e. the normal force at 
contact exactly balances the weight of the material con- 
tained in the half of the ruck for which S < L/2. We 
therefore find that a ruck sticks provided that a < ijlLq. 
This result is exact though, in general, the dependence 
of a on AL must be determined numerically. However, 
for small deformations, ([5]) may be used to relate a and 
AL and show that a ruck sticks only if 
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i.e. a ruck must be sufficiently large to stick. Qualita- 
tively speaking, this result remains valid even for large 
deformations, since the compressive force a is always a 
decreasing function of the imposed end-end compression, 
AL. Fig. [2] shows a regime diagram of the regions of 
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(/iLo, AL) parameter space for which a ruck sticks or 
slips. Fig. [2] also shows the results of experiments in 
which the critical value of AL at which rucks first stick 
was determined for different sheet lengths using a vari- 
ety of sheet materials and substrates. In each case, the 
value of £ g was measured experimentally using the loop 
test [22| and \x was measured from the angle of friction. 

We now move on to discuss the dynamic motion of 
a ruck. The experiments described here involved my- 
lar sheets (Goodfellow, Cambridge) of thicknesses h = 
125,250,350 and 500 and length 2 m. These sheets 
were laid horizontally on various different substrates (to 
investigate the role of substrate friction). The far end 
of the sheet is left free to move[23j and the near end at- 
tached to a vertical track. We generate a moving ruck 
by lifting the near end of the sheet vertically along the 
track (allowing the free end to slide in) and then moving 
the near end rapidly downwards. This leads to the for- 
mation of a localized ruck, which typically propagates at 
speeds ~ 1 m/s away from the near end. The size of the 
ruck may be controlled by lifting the near end to different 
initial heights. However, to render air resistance negligi- 
ble we ensure that the mass of air contained beneath the 
ruck is less than the mass of the ruck itself. This requires 
that d <C ph/p a where p a is the density of air. 

Propagating rucks were filmed using a high speed cam- 
era (Photron Fastcam) with a frame rate of 250 Hz. (See 
EPAPS Document No. [number will be inserted by pub- 
lisher] for a movie of a typical experiment.) The resulting 
images were analysed using ImageJ (NIH) to determine 
the ruck shape and the position of the peak of the ruck, 
Xp ea k, as functions of time. Typical experimental results 
are shown in fig. [3] The time dependence of X pea k, fig. Eh., 
shows that after some initial transient, the ruck moves at 
constant speed for a time before it slows down and then 
speeds up. This latter phase of the motion corresponds to 
the free end of the sheet beginning to slip relative to the 
substrate. The ruck profiles shown in fig. [3Jd show that 
while the ruck moves at constant speed its shape remains 
remarkably constant (though during the early transient 
and as the ruck slips away its shape does evolve). We 
focus on understanding this phase of the motion here 
and leave the early and late time behaviors to a future 
investigation. 

To derive the equations governing the dynamical mo- 
tion of a ruck, we use the approach of [HI, 0] albeit 
modified to include the vertical acceleration due to grav- 
ity. Resolving the stress within the sheet into horizontal, 
f( h ', and vertical, f^ v \ components, we may write the 
horizontal and vertical force balances as 

fi h) =X TT , f^=Y TT -l, (7) 

while the torque balance gives 

Oss - 4t 9 tt = f {h) Ys - f (v) X s , (8) 
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FIG. 3: (Color online) Experimental results for the dynamic 
propagation of a ruck in a mylar sheet (h — 125 pm). (a) 
The horizontal position of the peak of the ruck as a function 
of time (points). The speed in the steady state phase, c^/g£ g , 
is taken to be the gradient of the best fit line (dashed line) 
in the region where the shape is steady, (b) The shape of 
the travelling ruck at six instants of time. (The time of each 
profile is given by the position of the corresponding symbol 
in (a)). The solution of the heavy elastica equation @ with 
the same value of AL = 1.05 is shown by the solid curve. 

in which E is the Young's modulus of the sheet and T = 
t(g/£ g ) 1 ^ 2 is dimensionless time. 

The experimental results presented in fig. [3] suggest 
that there may be a travelling wave solution of eqns (JT]) 
and ijHJ). It is therefore natural to transfer into a frame 
moving with constant speed c. We introduce a new vari- 
able rj = S—cT, which enables us to integrate and de- 
termine the functions fM = F^(r?) and = F^(r]). 
Substituting these functions into {5| gives a single equa- 
tion for 9 = Q(i]): 

(1 - c 2 pg£ g /E) 9 W = - (<r + c 2 ) sin 9 + (L/2 - r?) cos 9. 

(9) 

In our experiments, the dimensional speed of the rucks, 
c 2 g£ g <C E/p, the speed of sound within the sheet. We 
may therefore neglect the difference between the prefac- 
tor of Q nv in © and unity so that © becomes exactly 
the heavy elastica equation with the eigenvalue a re- 
placed by an 'effective tension' a + c 2 . Thus, for a given 
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FIG. 4: (Color online) Main figure: Logarithmic plot of the 
effective tension, a + c 2 , (see text) measured experimentally 
for rucks in mylar sheets as a function of the dimensionless 
ruck height, 5 = d/£ g . Experimental points are shown for dif- 
ferent sheet thicknesses and substrates (hence different values 
of the friction coefficient fi) as shown in the legend. The val- 
ues measured of the friction coefficients are [i — 0.13 (copper), 
fj, = 0.41 (carpet) and /j, — 0.16 (wood). The dependence of 
a + c 2 on 5 determined from the numerical solution of @ is 
plotted as the solid curve while the asymptotic result (|10l) is 
plotted as the dashed line. Inset: Linear plot of dimensional 
raw data showing the dependence of ruck speed c^J g£ g on 
ruck height d. 



value of AL, the shape of a steadily moving dynamic 
ruck must be exactly that of the static ruck with the 
same value of AL. This point is illustrated in fig. [3)3 
where we see that while the ruck is propagating with a 
constant speed, its shape is indistinguishable from that 
of a static ruck with the same value of AL. 

It is a simple matter to solve numerically the eigen- 
problem |9]) with boundary conditions analogous to f4|. 
The linearized problem may be solved analytically giving 
the effective tension in terms of the ruck height S = d/£ g 



3.962<T 1/2 . 



(10) 



The presence of the combination a + c 2 in ([9]) and 
(fl~0| indicates that the ruck speed c is not uniquely de- 
termined by the solution of the eigenproblem ([9]) with 
boundary conditions ([4]). This is demonstrated by the 
inset of fig. |3J which shows raw results for the speed, 
Cy/g£ g , as a function of ruck height, d. However, we note 
that the slipping away of the ruck enables us to estimate 
the value of a — at the commencement of slipping, a 
must exactly balance the maximum static friction force 
that can be generated by the material remaining between 
the ruck and the free end. Fig. [5] shows non-dimensional 
results for er+c 2 as a function of ruck height S determined 
using this experimental procedure. Also plotted in fig. [1] 
are the theoretical predictions obtained by solving the 
full problem numerically (solid curve) and the result of 
the linear analysis (fTT))) (dashed line). These show that 



we obtain good quantitative agreement between theory 
and experiment, though inaccuracies in determining the 
onset of sliding limit this agreement. 

In this Letter we have considered the properties of 
static rug rucks and shown that friction allows sufficiently 
large rucks to remain once the initial compression is re- 
moved. We have also considered the inertial dynamics 
of rucks, complementing a recent study of the creeping 
motion of a ruck on an inclined plane [2l| . A result of par- 
ticular interest is that large rucks generally move more 
slowly than smaller ones. If colliding rucks aggregate 
this fact would drive a population of travelling rucks to 
form one large, slow moving ruck. Future work will fo- 
cus on determining whether this is qualitatively the same 
in populations of slip pulses in geophysical settings and 
could in turn lead to alternative mechanical rationaliza- 
tions [IHl of the statistics of earthquakes. 

We thank L. Mahadevan for bringing a related problem 
2l| to our attention in 2003, arousing our interest in 
rucks. D.V. is supported by the 1851 Royal Commission. 
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